iV-conserving Bogoliubov vacuum of a two component Bose-Einstein condensate: 
Density fluctuations close to a phase separation condition 



Bartlomiej Oles and Krzysztof Sacha 
Marian Smoluchowski Institute of Physics and Mark Kac Complex Systems Research Center, 
Jagiellonian University, Reymonta 4, 30-059 Krakow, Poland 
(Dated: March 26, 2008) 

Two component Bose-Einstein condensates are considered within a number conserving version of 
the Bogoliubov theory. We show that the Bogoliubov vacuum state can be obtained in the particle 
representation in a simple form. We predict considerable density fluctuations in finite systems close 
to the phase separation regime. We analyze homogeneous condensates and condensates in a double 
well potential. 
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I. INTRODUCTION 

Bose-Einstein condensate (BEC) is a unique state of 
a many particle system where, ideally, all particles oc- 
cupy the same single particle state. It is obviously pos- 
sible for bosons only, and experimentally it can be real- 
ized in ultra-cold dilute atomic gases Since the first 
experimental realization numerous different phenomena 
involving BEC have been investigated and nowadays it 
is also possible to obtain mixtures of BECs or even mix- 
tures of ultra-cold bosonic and fermionic gasses Q ■ Two- 
component BEC [H can reveal number of interesting phe- 
nomena, e.g., phase separation [1,3, self-localization Q, 
condensate entanglement [(| @] or internal Josephson ef- 
fects 0. 

In an infinite homogeneous system the phase separa- 
tion occurs abruptly once interactions reach their critical 
values 0, [|| . In the present paper we show effects of a 
finite system. That is, in a finite box there is a region 
close to critical values of the coupling parameters where 
substantial density fluctuations can be observed. 

A standard theoretical description of a single conden- 
sate and condensate mixtures starts with the mean field 
Gross Pitaevskii equations @ that provide estimates for 
ground states and collective excitations of a system but 
under an assumption that it is described by perfect con- 
densate product states. Particle interactions, however, 
can lead to substantial depletion of the condensates [l| 
and in order to obtain a more realistic picture usually 
a Bogoliubov theory is applied, which allows one to de- 
scribe small quantum corrections to the mean field solu- 
tion P, [13, ELIpl- The key idea of the original Bogoli- 
ubov theory [T(| (usually used in the BEC field) is the 
U(l) symmetry breaking approach where the atomic field 
operator is assumed to have a nonzero expectation value. 
This coherent state necessarily involves superposition of 
states with different numbers of atoms, an assumption 
very far from experimental reality. Moreover, careful 
analysis of the original theory shows that the Bogoliubov- 
de Gennes equations correspond to an eigenvalue problem 
of an operator which is not diagonalizable and the theory 
must break down after a finite time [ll|, EH • 

To overcome these drawbacks we employ a number 
conserving version of the Bogoliubov theory, which has 



been presented by Castin and Dum [ll| (see also [12 ]) 
for a one-component BEC and generalized to a two- 
component system by S0rensen pj, and analyze the Bo- 
goliubov vacuum state of a mixture of two BECs. The 
two theories should give the same physical predictions for 
large particle numbers. There are, however, examples of 
systems where the TV-conserving theory works in a regime 
of the standard theory breakdown . The Bogoliubov 
vacuum is usually obtained in the quasi-particle repre- 
sentation where quantum depletion, i.e. the number of 
particles occupying non-condensate modes can be easily 
calculated for a given system p], QjJ |Tl|, [l2j| ■ To gain in- 
sight into the form of the ground state of the system, we 
derive the Bogoliubov ground state in the particle rep- 
resentation. This enables us to perform simulations of 
density measurements in single experiments [HI, [lj| [l7| ■ 



For a single condensate the Leggett ansatz of the vac- 
uum for translationally invariant systems [l8| has been 
shown to be valid in any inhomogeneous condensates in 
[Til . flU ] . In the present paper we show that the ansatz 
can be used also in the two-component case even in the 
presence of the inter-species interaction. The obtained 
Bogoliubov vacuum state is then used in an analysis of 
density fluctuations in 3D homogeneous condensates and 
in condensates trapped in a double well potential. It 
turns out that vicinity of the critical point for the phase 
separation is especially interesting because the fluctua- 
tions there become considerable. 



The paper is organized as follows. In Sec. II we present 
the solution for the Bogoliubov vacuum state in the parti- 
cle representation, derived within the number conserving 
version of the Bogoliubov theory. In Sec. Ill we describe 
a procedure used later to perform density measurement 
simulations. The theory is applied to the analysis of ho- 
mogeneous condensates in Sec. IV and to the double well 
problem in Sec. V. We conclude in Sec. VI. Short reminder 
of the Bogoliubov theory Q is presented in the Appendix 
A and details of the derivation of the Bogoliubov vacuum 
state in the particle representation are presented in the 
Appendix B. 
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II. BOGOLIUBOV VACUUM STATE 

We consider a two component Bose-Einstein conden- 
sate formed by a mixture of two kinds of atoms (or the 
same atoms in two different internal states), i.e. N a 
atoms of type a and N b atoms of type b. The Hamil- 
tonian of the system reads 
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2m b 
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where m a , mj arc particle masses, V a (r), V b (r) stand for 
the trapping potentials and 
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where a a , a bl a ab are the scattering lengths. The number 
conserving Bogoliubov theory 0, \1M assumes the follow- 
ing decomposition of the bosonic field operators 



tp b (r) = <j>oo(r)b + Sipb(r), 



(3) 



where we separate the operators do and bo that anni- 
hilate atoms in modes 4> a o and <p b o, respectively, which 
are macroscopically occupied by atoms. That is, for the 
states we are after 



N a , 



N b 



(4) 



Corrections 8ip a and Sip b are thus supposed to be small 
and we may perform expansion of the Hamiltonian in 
powers of 5ip a and Stp b . In the zero order, condition for 
energy extremum in the 4>ao an d </>&o space leads to cou- 
pled Gross-Pitaevskii equations 



where 
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Ma, 



Mb, 
(6) 



(with chemical potentials fi a and n b ) that allow us to find 
single particle modes macroscopically occupied by atoms. 
The first order terms of the Hamiltonian disappear. In 
the second order one obtains an effective Hamiltonian 



which, employing the Bogoliubov transformation, can be 
written in a diagonal form 



(7) 



n£"+" 



where the sum goes over the so-called family "+" solu- 
tion of the Bogoliubov equations (see Appendix A). The 
quasi-particlc annihilation operators arc defined as: 



«|A )-«|At) + (^|A 6 )-(^|Ah, (8) 



where 
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The wavefunctions are solutions of the 

Bogoliubov equations corresponding to eigenvalue E n 
(see Appendix A). Let us now switch to our results. 

The Bogoliubov vacuum state |0s) is an eigenstate of 
the effective Hamiltonian that is annihilated by all quasi- 
particlc annihilation operators, 



Cn\0 B ) = 0. 



(10) 



Other eigenstates can be generated by acting with the 
quasi-particle creation operators & n on the Bogoliubov 
vacuum. The quasi-particle representation is thus natu- 
ral to represent the system eigenstates within the Bogoli- 
ubov theory. It is also suitable to obtain low order cor- 
relation functions. However, to get predictions for den- 
sity measurements, i.e. to simulate measurements of all 
atom positions, the particle representation turns out to 
be much more convenient. 

In the Appendix B we show that the Bogoliubov vac- 
uum state can be written in the particle representation 
in the following simple from 
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where 
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The particle creation operators (6^) create particles in 
modes (\> aa {4>ba ) that are eigenstates of the single particle 
density matrices, and dN® {dN^) are the corresponding 
eigenvalues, i.e. 

(0 B \^ a (r)Mr')\0 B ) a A a C o (r)0 aO (r") 
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and similarly for (ObI^M^o^OIO-b)- The operators f a 
and Tb are denned as: 
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The presented solution (fTT]) is self-consistent provided the 
Ta,b operators are diagonal in the basis of the eigenvectors 
of the single particle density matrices. In the following 
sections we show examples of a spatially homogeneous 
system and BECs in a double well potential, where this 
indeed is the case. 



III. DENSITY MEASUREMENT 

Average particle density corresponds to the reduced 
single particle density which can be easily calculated 
within the Bogoliubov theory The average density 
means an averaged picture obtained by collecting out- 
comes of the density measurement in many experimen- 
tal realizations of a system in the same quantum state. 
Even at zero temperature a many body system, can re- 
veal density fluctuations and a single photo of the system 
may be sig nificantly different from the averaged picture 

In order to perform density measurement simulations 
we generally need a full many body probability den- 
sity. As the number of particles grows, however, using 
this density quickly becomes a very formidable task. In- 
stead one may use a sequential method proposed by Ja- 
vanainen and Yoo (l~5| . It relies on a choice of a position 
of a subsequent atom with the help of a conditional den- 
sity probability which takes it into account that previ- 
ous atoms have already been found at certain positions. 
Note that, since this method requires acting with particle 
annihilation operators on the Bogoliubov vacuum, using 
the Bogoliubov state in the quasi-particle representation 
would require inversion of the nonlinear transformation 
©. Having the state (fTTj) we avoid this problem. 

In practice the sequential method fl5| can be used if 
only one (or few) non-condensate modes are important. 
If many modes are relevant we should e.g. switch to an 
approximate method [T(| . Suppose there are M a and Mb 
modes where 
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Then results of the density measurements corresponding 
to a state of the form (fTTj) can be approximated by [l6[ 
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where 



Vba (r) 
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(17) 



and real parameters q aa and qb a have to be chosen ran- 
domly, for each experimental realization, according to a 
Gaussian probability density 
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The replacement (fT7|) is essential because it makes all 
eigenvalues of the T a ^b operators non-negative and allows 
writing the Bogoliubov vacuum state of the form (fTTj) 
as a Gaussian superposition over condensates which, in 
turn, leads directly to the predictions (|16[) (l6| . 



IV. HOMOGENEOUS CONDENSATES 

The two component homogeneous condensate is an ex- 
ample of a Bose system where the Bogoliubov theory 
gives analytical results even in the presence of a pro- 
cess which transfers atoms between the two components 
(Rabi or Josephson coupling In numerous papers the 
quasi-particle excitation spectrum is analyzed as well as 
its dynamical instability leading to the phase separation 
[J, Qj| . In the present publication we fix the number of 
atoms in each component and study the Bogoliubov vac- 
uum state in the particle representation for interaction 
parameters approaching the phase separation condition. 

Suppose we deal with a condensate mixture in a box 
of L x L x L size with periodic boundary conditions and 
all interactions are of repulsive character, i.e. g a , gb, 
g > 0. The ground state solution of the Gross-Pitaevskii 
equations reveals condensate wavefunctions 



1 



»a0 



'Ml 



and chemical potentials 



Ma = gaPa + gPb, Mb = 9bPb + gPa, 



(19) 



(20) 



where p a ^ = N a ^/L 3 are densities of a and b compo- 
nents. For a homogeneous system it is appropriate to 
switch to the momentum space and look for the solution 
of the Bogoliubov-de Gennes equation in the form 
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Then, one obtains two quasi-particles for each k with 
energies 0,i, 
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where 
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and the operators a! 
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create atoms in 



the modes (|28]>. 

In the case of the infinite box (i.e. for L — > oo) if _g 2 > 
gagb uniform solutions of the Gross-Pitaevskii equations 
become unstable — mixing of the a and b components 
is not energetically favorable and the phase separation 
occurs [H H, Ojjj]. It manifests itself in the appearance 
of an imaginary eigenvalue in the Bogoliubov spectrum 
(f2"2"| . In the case of a finite box the minimal value of 
the momentum becomes 2tt/L and the condition for the 
phase separation is modified, 
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Note that in the finite box the momenta are discrete 

+ n y e y + n z e z ), (27) 



k = — (n 7 



where n x ,n y ,n z are non-zero integers. 

The reduced single particle density matrices are diag- 



which shows that for a finite system the minimal value 
of the parameter g leading to the phase separation has 
to be greater than the corresponding value for L — > oo. 

In Fig. Q] we show average numbers of atoms depleted 
to the modes (|28p and values of the corresponding X k ' b , 
Eqs. (f3"0"|) , far from the phase separation. The data cor- 
(26) res P orL( i to a mixture of 87 Rb atoms in two different in- 
ternal states, N a = 5000, N b = 20000, L = 50 pm, 
a a = 108. 8ao, a b = 109. lao and a a b = lOaoi where ao 
is the Bohr radius [20]. The value of the latter scat- 
tering length can be adjusted by means of a Feshbach 
resonance Q. To obtain predictions for atomic density 
measurements one has to change phases of the modes, see 
(fP7]) . which in the present case of all negative A^' fe leads 



onal in the e /vL 3 basis. However, in order to have 
the F a operators also diagonal we have to switch to the 
basis 
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and similarly for the <p a j: c and f bkc modes. Due to the 
fact that </) a o and (j) b o are real and all the ip modes are 
purely imaginary we obtain, see (|16p . 



Then the Bogoliubov vacuum state in the particle repre- 
sentation reads 
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A b j2 and A^ fc « N a ^ the 
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density fluctuations turn out to be negligible, i.c 
density remains almost perfectly flat. 
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FIG. 1: Panel (a) presents condensate depletion, i.e. number 
of atoms depleted from a condensate wavefunction, dN£' = 
2 + ( v k',-) 2 y to modes (JgHJ). Panel (b) shows the corre- 
sponding values of the parameters A^ ,!) , Eqs. (|30p . Circles are 
related to the BEC component a while crosses to the com- 
ponent b. The results correspond to the parameters far away 
from the phase separation, i.e. N a = 5000, Nt, — 20000, 
L = 50 fim, a a — 108. 8ao, a,b = 109. lao and a a b = lO.Oao, 
where ao is the Bohr radius. 



Figure [2] shows similar data as Fig. Q] but for a a b = 
193. 9a , i.e. chosen so that the phase separation condi- 
tion in an infinite system would be already fulfilled but 
it is still not fulfilled in the case of the finite box. The 
numbers of atoms depleted are not dramatically greater 
than the ones in the case considered previously but now 
some values of A^' h become positive. The latter has dra- 
matic consequences for density fluctuations because the 
modes ip corresponding to the positive X^' b are real and 
their contributions to the atomic density are of order 

of \J Aa' h /N a fi. Indeed, the predictions for atomic den- 
sity measurements (neglecting contributions of order of 
Kl b l N a,b) show that 



FIG. 2: The same as Fig. [T] but for conditions close to the 
phase separation, i.e. a a t = 193. 9ao. 
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FIG. 3: Solid curves show predictions for results of den- 
sity measurements of atoms belonging to the a component 
- the densities integrated over y and z directions, i.e. 
L J a a (r) dydz, are shown. Thick solid line denotes the den- 
sity averaged over many experimental realizations (i.e. the re- 
duced single particle density) while dashed lines indicate the 
average density plus/minus standard deviation that is equal 
; to 2^/Al/Na fa 0.1 (where k = 2n/L). 
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where y"V runs over modes corresponding to positive 
A^' b only. In each experimental realization one has to 
choose g ojEaI g a£c , q bJ : s and g fejSc randomly according to 
the probability density ([T5f . In Fig. Owe show a few 
examples of the simulations for atoms belonging to the a 
component together with the averaged result (which cor- 
responds to the reduced single particle density) — the 
figure presents the densities integrated over y and z di- 
rections. Despite the small number of atoms depleted 
to the lowest momentum mode (~ 0.3%) the changes of 
the density are of order of 10%. Standard deviations of 
the largest scale density fluctuations (i.e. corresponding 
to the quasi-particles with the momentum k = 2-k/L) 
behave like 
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where a ab is the critical value for the phase separation. 
Exactly at the critical point the Bogoliubov theory breaks 
down, which is indicated e.g. by the divergence of the 
fluctuations in Eq. ([35]). For the parameters chosen in 
Fig. [2][3] we are, however, sufficiently far away from the 
critical point so that the predictions on the basis of the 
Bogoliubov theory are reliable. 

From the experimental point of view it is important 
that the density fluctuates on a scale of the order of L. 
One may use low resolution in the density measurements 
so that statistical fluctuations will be practically elimi- 
nated and the only density modulations will correspond 
to the fluctuations considered here. 

In the example considered, the range of a a b where one 
deals with positive A?' 6 is about 10 x Bohr radius, which 
should be wide enough to enable experiments with the 
density fluctuations ([331) . 

Note that the structure of the Bogoliubov vacuum 
state ([2T)|) shows that there are no correlations between 
atoms belonging to the different components. Recently 
in Ref. [2l|, Bogoliubov vacuum of the form of (jTTj) has 
been used in a two component system in the case when 
the inter-species interaction is absent, g — 0, and the 
components become fully independent. Our analysis in- 
dicates that the Bogoliubov vacuum possesses the same 
form even in the presence of the inter-species interaction. 
Of course, the interaction changes the Bogoliubov modes 
and influences the values of A2/ fc . 



V. DOUBLE WELL 

In the present section we will consider a simple model 
where there are analytical solutions within the Bogoli- 
ubov theory both in the miscible and in the phase sepa- 
rated regime. 



Let us consider a two component Bosc-Einstein con- 
densate in a one-dimensional symmetric double well po- 
tential under an assumption that the Hilbert space of the 
system is restricted to ground states in each well only (i.e. 
within the two mode approximation). For experimen- 
tal realizations of the double well problem see [22|]. The 
Hamiltonian of the system, if we choose real functions as 
the ground states in the two wells, reads 
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where the A\ [B\) operator annihilates an atom belong- 
ing to the component a (b) in the first well and the A 2 
(B 2 ) operator annihilates an atom of the a (b) component 
in the other well. For the calculations we have chosen po- 
tential wells situated at x = — 2, 2 and with such widths 
that the ground states of the wells arc 
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The parameter Q stands for the frequency of the tun- 
neling of atoms between the two wells and U and U a b 
describe intra- and inter-condensate interactions. We 
will consider the case of symmetric BEC components, 
i.e. N = N a = N b = 1000, Q = 5000, U = 1 and U ab > 
(i.e. all interactions are of repulsive nature), but similar 
analysis can be easily performed for N a ^ N b and for tun- 
neling frequencies and intra-species interactions different 
for both components. 

We would like to mention that in the limit of large at- 
tractive interactions a two-component entanglement has 
been found in the system [23| . 



A. Mean field solutions 

For U ab smaller than the critical value 

Ulb = Jf+ U > 



(38) 



the ground state solution of the Gross-Pitaevskii equa- 
tions (J5[) reveals both condensates symmetrically located 
in the double well potential: 

<l>a,o(x) = <l>bo(x) = -j= [ipi(x) + ip 2 (x)] , (39) 

i.e. we are in the miscible regime. However, if the pa- 
rameters of the system fulfil U ab > U^ b the solution ([3"9"[) 
is unstable - the spatial overlap of the atomic clouds of 



7 



the different components becomes energetically not fa- 
vorable (see Fig. 2^) and the phase separation begins. 
The condensate wavefunctions arc then: 
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Note that in the phase separation regime there are two 
ground state solutions, for exchanging a <-> /3 in (|40[) one 
obtains another solution of the Gross-Pitaevskii equa- 
tions. On the basis of these two solutions two different 
Bogoliubov vacuum states can be obtained. In the fol- 
lowing we will show that sufficiently far away from the 
critical point the ground state of the system can be ap- 
proximated by preparing a superposition of the two Bo- 
goliubov vacuum states. 



B. The Bogoliubov vacuum — miscible regime 

The solution of the Bogoliubov-de Gennes equations 
reveals two quasi-particles corresponding to energies: 
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and modes: 
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where 



and 



4>ai{x) = 06i (x) = -j= [ipi(x) - i>i{x)} 
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is the normalization factor. 

The Bogoliubov vacuum state can be written in the 
form (fTTj) using 
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FIG. 4: Panel (a): dashed line shows mean field energy of 
the symmetric Gross-Pitaevskii solution (|39[1 while solid line 
denotes the energy of the asymmetric solution (|40[) which ap- 
pears in the phase separation region. Panel (b) shows de- 
pletion of the condensates while panel (c) the corresponding 
values of X a ' b . Note that, since we consider symmetric in- 
teractions, the depletions are equal for both components and 
A Q = X b . 



The parameters (|47[) and the depletion arc shown versus 
U ab in Fig. m 

As can be seen from the quasiparticle excitation en- 
ergy (|43|) . when the condition U a b > U^ b is fulfilled, the 
spectrum is no longer real and the solution (|39[) becomes 
unstable. The ground state of the Gross-Pitaevskii equa- 
tions shows separation of the atomic clouds of the a and 
b components. 



C. The Bogoliubov vacuum — phase separation 
regime 



Now the quasi-particle excitation energies are 
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1/2 



± 2W ab a(3 



(48) 
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and the quasiparticle modes, corresponding to the con- 
densate wavefunctions (|40[) . are proportional to 



Pbi(x) 



(49) 



We skip here rather long expressions for the quasiparti- 
cles since they can be easily obtained with the help of the 
Bogoliubov transformation. Behaviour of the component 
a (or b) depletion and of the parameters A°, X b versus U ab 
is depicted in Fig. |U 

A Hamiltonian of a BEC system in a symmetric double 
well potential is invariant under the parity inversion, i.e. 
if we reverse coordinates of all particles the Hamiltonian 
does not change. It implies that eigenstates of our system 
(unless there is a degeneracy) must be also eigenstates of 
the parity operator. In the phase separation regime the 
Gross-Pitaevskii solutions (|40|) are neither even nor odd 
functions and the corresponding Bogoliubov vacuum, 



N/2 
,-, N/2 



(50) 



is not an eigenstate of the parity operator. Exchanging 
a with (3 one obtains another Gross-Pitaevskii solution 
and another Bogoliubov vacuum state. A proper parity 
state can be obtained by preparing a superposition of the 
two states 



|0 B ) + |0fl(a «->£)>. 



(51) 



The state ([ST]) is a good approximation for the ground 
state of the system if we are not very close to the critical 
point. That is, if we increase U ab for fixed N, the states 
|0b) and |0s(a <-> 0)) very quickly become practically or- 
thogonal, and the sooner it takes place, the greater N we 
choose. Note that even in the regime of these states be- 
ing orthogonal, the corresponding mean field states (|40|) 
need not be orthogonal at all. If we are, however, far 
away from the critical point also the Gross-Pitaevskii 
solutions have zero overlap, i.e. {4> a o\4>bo) ~ 0. Then 
the state (f5"T|) is a a Schrodinger cat state [24| which is 
strongly vulnerable to atomic losses — loss of a small 
number of atoms is sufficient to distort completely the 
coherent superposition in (|5ip . 

At the critical point the Bogoliubov theory breaks 
down because higher order terms become dominant. 



D. Density fluctuations 

We see in Fig. |4] that approaching the critical point 
in the miscible regime the depletion of the condensates 
and the Bogoliubov vacuum parameters A a ^ grow. Very 
close to the critical point the depletion is very large and 
the Bogoliubov theory can not be longer applied. How- 
ever, in the vicinity of the point there is a regime where 



the depictions arc very small compared to the total par- 
ticle numbers and the parameters A a ^) are positive. The 
modes (f) a Q , (f>bo , <p a i and 4>b\ are real and the appearance 
of the positive A a ' b indicates (similarly as in the previous 
section) that density fluctuations become considerable, 
i.e. of order of ^k a > b /N, 

a a(x) ~ 4> 2 a0 {x) + 2-?=<j) aQ (x)(f) al (x), 
v TV 



a b (x) ~ 0to(x) + 2—=4> m (x)<j) bl (x), 

v Jy 



(52) 

where q a and q b have to be chosen randomly according to 
(|18p in order to get predictions for the results of the den- 
sity measurement in different experimental realizations. 

A few examples of the density measurements of the a 
component atoms in the miscible regime are shown in 
Fig. [5] for U^ b — U a b = 0.005. Standard deviations of the 
density fluctuations behave like: 



(53) 



(u c ab - U ab ) 



1/4 1 



i.e. similarly as in the homogeneous case, see (|35[) . 

On the other side of the critical point, i.e. in the phase 
separation regime, we deal with a state of the form (|5Tj) 
which for Ul b -U ab = -0.005 and N = 5000 is a good ap- 
proximation for the ground state of the system, indeed 
I (0b|0b(q! «-> (3)) | 2 is of order of 10 -8 . One may also 
expect (similarly as in the miscible regime) substantial 
density fluctuations. However, for a superposition of the 
Bogoliubov vacuum states (|5ip we cannot simulate den- 
sity measurements with the help of the method described 
in Sec. III. 



VI. CONCLUSIONS 

We have considered a number conserving version of the 
Bogoliubov theory for a two component Bose-Einstein 
condensate, with the fixed number of atoms in each com- 
ponent. We have shown that the Bogoliubov vacuum 
state can be written in the particle representation in a 
simple form, provided that eigenstates of the reduced sin- 
gle particle density matrices diagonalize the operators 
(|14[) . Having the Bogoliubov vacuum in the particle rep- 
resentation one can easily obtain predictions for density 
measurements in single experiments. 

The introduced formalism has been applied to the anal- 
ysis of a two component homogeneous condensate and a 
two component condensate in a double well potential. In 
finite homogeneous systems, when parameters of the sys- 
tem approach a phase separation condition, considerable 
density fluctuations appear before the system becomes 
unstable. This behaviour is different than in infinite sys- 
tems, where the phase separation happens abruptly. The 
range of the parameter values where the substantial fluc- 
tuations are observed indicates that the results presented 
here can be verified experimentally 
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0.5 




FIG. 5: Density fluctuations of the a component in the mis- 
cible regime close to the critical point, i.e. for U a b = 5.995. 
Solid line corresponds to the mean field solution, dashed and 
dotted lines are two examples of possible realizations of the 
experiment. The standard deviation of the density fluctuation 
is 

In the case of condensates in a double well potential 
we are able to describe the system in a vicinity of the 
critical point both in the miscible condensates regime 
and in the phase separation region. Considerable density 
fluctuations can be expected if the parameters approach 
the critical values. 
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APPENDIX A 



We begin with a short reminder of the results of the 
number conserving version of the Bogoliubov theory. Fol- 
lowing 0] we will perform the perturbation expansion of 
the hamiltonian. The decomposition ([3]) allows us to ex- 
pand the Hamiltonian in powers of small operators 8?p a 
and 5ipb- As mentioned in Sec. HH minimizing the en- 
ergy of the system in the zero order we obtain coupled 
Gross-Pitaevskii equations j5]) that allow us to find the 
condensate wavefunctions a o and (f> b Q. The first order 
terms of the Hamiltonian disappear and in the second 
order we obtain an effective Hamiltonian 




(A-l) 



where 



gV 'N a N b Q a 4>aQ(t)wQl \ 

-gVWbQ* a 4>*ao4>boQ* b 
g b N b Q b 4 Q* b 

-H b GP -g b N b Qt\M 2 Ql I 

Note that action of the A Q and Af, operators preserves 
numbers of atoms in the system. Diagonalization of 
the effective Hamiltonian amounts to solving the eigen- 
equation for the non-hermitian operator C (i.e. the 
Bogoliubov-de Gennes equations). 

The C operator possesses two symmetries (similarly 
to the symmetries of the original Bogoliubov-de Gennes 
equations [ill]). 

u\Lu\ = —£*, 



gaN a Q a 4> 2 a0 Q* a 
-H^ p -g a N a Q*M 2 Q* a 

g\J N a N b Qb(j>ao4>bC>Q*a 
-g \/N a NbQ* b (f>a0 060 Q*a 



g\/ NaNbQa^aO^boQb 

-gVNMQtKoKoQb 

H b GP +g b N b Q b \<Pbo\ 2 Qb 
~g b N b Q* b cj)l 2 Qb 
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where 



"1 



and 





u 3 Cu 3 


CTl 


\r)- 


U 


i / 


(! 


!)■ 







"3 = 



<?3 








CT3 



(A-5) 



(A-6) 



CT3 



1 

-1 



are the first and third Pauli matrices, respectively. Sup- 
pose that all eigenvalues of the C operator are real. The 
symmetries (|A-5[) imply that if 



( K> \ 



(A- 



is a right eigenvector of the C with eigenvalue E n , then 
|^) = U3|^^) is a left eigenvector of the same eigenvalue 



E n . and ui^f*) is a right eigenvector with eigenvalue 
—-En- 
There are four eigenvectors of C corresponding to a 
zero eigenvalue, 



(A-9) 



The other eigenstates of the C operator we divide into 
two families "+" and "— ", according to 





( 10a) \ 


( ^ 




( o ^ 




( o \ 







m 












(A-7) 
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^ o J 


V o j 




^ o j 




^ \<t>l) / 



(A-10) 



Having the complete set of the eigenvectors of the C we 
obtain an important completeness relation 



/ l<) \ 



E 



/ K*) \ 



V M) ) 

( |0ao)(0ao| 

ICoXCol 

o o 



;0 + E 

n£"+" 



(-(CI, «*!.-(«» I, <«n I) 



v 



o 











(A-ll) 



The eigenvectors of the L operator define the Bogoliubov 
transformation 



(K\ 

A 
A b 



E 



V < J 



E 



( <* \ 



(A-12) 

3|) fulfill the bosonic 
commutation relation [c„,cL] ~ 5 n . n '- Employing the 
Bogoliubov transformation we obtain the effective Hamil- 
tonian in a diagonal form l[7]). 



where the quasi-particlc operators 



APPENDIX B 

The Bogoliubov vacuum state \0b) is an eigenstate of 
the effective Hamiltonian ([7]) that is annihilated by all 



quasi-particlc annihilation operators. Let us show that 
the Bogoliubov vacuum can be obtained from the particle 
vacuum by applying some particle creation operators cfi a 
and dt 



10, 



dl) (dt 



10), 



(B-l) 



where we require that d\ and dl commute with all quasi- 
particlc annihilation operators [161 ]. 



[CnJt] 



= o, 

= 0. 



(B-2) 



Then the state (|B-1[) is indeed annihilated by all quasi- 
particle annihilation operators, 



(4) 



Ma / \ Mi 



|0) 



. v M„ / N Mi 



6 y 2»|0) = 0. 

(B-3) 
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The set of equations (|B-2[) is solved by the particle 
creation operators in the form [l6|, [l8[ 



ne"+" 



OR* I 



(B-9) 



a.. 0=1 

oo 



4 = bihi+ J2 zipVJl* 



a, 0=1 



(B-4) 



where 1" and lj_ are the identity operators in the sub- 
spaces orthogonal to the condensate wavefunctions <f> a o 
and 4>bo, respectively. Comparing Eq. (|B-8[) with 



where dj, (b\) are bosonic particle creation operators that 
create atoms in modes 4> aa (<fib a ) orthogonal to the con- 
densate wavefunction (j) a0 (fibo)- %a0 and Z h a(i are sym- 
metric matrices to be found. 

Substituting the ansatz (|B-4[) into (|B-2[) we obtain 
equations: 



(0 B \4,t(r)Mr')\0 B ) = N a 4> aQ (r)<t> M (r') 



0=1 

OO 

i<\tia) = £<^|<M4- (B-5) 
0=1 

which, when multiplied by ((/> 07 |w°) and ((fib-yln 1 ^) , respec- 
tively, and summed over n, are transformed to 



■ £ <(fK*(f'!jB-io) 

ne"+" 



we see that the U a operator is a sum of a part of the 
reduced single particle density operator corresponding to 
the subspace orthogonal to the condensate wavefunction 
(j) a o and the identity operator 1". Similar statement is 
true in the case of the b component. Thus, if we choose 
as a basis (j) a a (</>b a ), the eigenstates of the single particle 
density matrix, we get the U a (Ub) operator in a diagonal 
form. Then one obtains immediately the solutions for the 
Z a Ji matrices, i.e. 



(0a 7 |fa|0aa> = zl (0a 7 | U a 1 a/3 ) Z% a 
0=1 

OO 

(^\T b \r ba ) = J2(<t»n\Ub\MZ b p a , (B-6) 
0=i 



where 



u a = J2 K>«i» u » = J2 \OM\> ( B - ? ) 

n£"+" n£"+" 

and T a: fc arc defined in (H3J). The completeness relation 
(|A-1 II) implies that the T a and Tb operators are symmet- 
ric and that 



(B- 



rya 



|r«| 



'a0l 



dN% + l 



7 b 



dN b „ + 1 



> (B-ll) 



where dN£ ,b are the eigenvalues of the single particle den- 
sity matrices, that is numbers of atoms depleted from 
the condensate wavefunctions to other eigenmodes. The 
Z a/3' zb oL0i (<Paa\Ta\4>*af3) and (<fibg |T b | fe) m atrices are 
symmetric. Thus, in order the ansatz (|B-4|) to be self- 
consistent the r a (, operators have to be also diagonal in 
the basis of the eigenvectors of the single particle density 
matrices. In Sec. IIVI and IVl we show examples where in- 
deed this is the case. Final form of the solution for the 
Bogoliubov vacuum state in the particle representation 
is presented in (jTTJ) . 
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